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ABSTRACT 



Context. The detailed comparison between observations and simulations of ram pressure stripped spiral galaxies in the Virgo 
cluster has led to a three dimensional view of the galaxy orbits within the hot intracluster médium. The 3D velocities and Mach 
numbers derived from simulations can be used to derive simple Mach cone geometries for Virgo spiral galaxies. 
' Aims. We search for indications of hot gas within Mach cones in X-ray observations of selected Virgo Cluster spiral galaxies 
CLil (NGC 4569, NGC 4388, and NGC 4501). 

^ Methods. Low-resolution maps of diffuse extended emission and X-ray spectra from XMM-Newton observations are presented. 
^ Gas donsities and temperatures are derived from the X-ray spectra. 

Results. We find extraplanar diffuse X-ray emission in all galaxies. Based on the 3D velocity vectors from dynamical modelling a 
■ simple Mach cone is fitted to the triangular shape of NGC 4569's diffuse X-ray emission. Assuming that all extraplanar diffuse 
I ' ' X-ray emission has to be located inside the Mach cone, we also fit Mach cones to NGC 4388's and NGC 4501's extraplanar 
, X-ray emission. For NGC 4569 it is hard to reconcile the derived Mach cone opening angle with a Mach number based on the 
' ^ ' sound speed alone. Instead, a Mach number involving the Alfvénic speed seems to be more appropriate, yielding a magnetic 
\ field strength of ~ 3-6 ¡iG for a intracluster médium density of n ~ lO^** cm""^. Whereas the temperature of the hot component 
I . of NGC 4569's X-ray halo (0.5 keV) is at the high end but typical for a galactic outflow, the temperature of the hot gas tails 
^ ■ of NGC 4388 and NGC 4501 are significantly hotter (0.7-0.9 keV). 

, Conclusions. In NGC 4569 we find direct evidence for a Mach cone which is fiUed with hot gas from a galactic superwind. We 
suggest that the high gas temperatures in the X-ray tails of NGC 4388 and NGC 4501 are due to the mixing of the stripped 
ISM into the hot intracluster médium of the Virgo cluster. 
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1. Introduction galaxies with long Hi tails in intermediate to lew galaxy- 

^ ' density regions (0.6-1 Mpc in projection from M87). The 

There is ampie evidence that the cluster environment ^^-^^ pointing roughly away from M87. Assuming 

modifies the gas distribution of spnal galaxies. Especially ^^¿-^^ g^j^^y ^^^-^^ ^^-^ suggests that these tails are due to 

m the Vngo cluster, where the gas distribution of spiral pressure stripping. Therefore, ram pressure stripping 

galaxies can be observed m detail by current mterferomet- ^j^^^^y ^^^^^^ ^^^^^ gpi^^l galaxies around the cluster 

ric observations, radially truncated gas disk are frequently yirial radius. The further evolution of a galaxy depends 

observed (Cayatte et al. 1990, Chung et al. 2009). The critically on its orbit (see, e.g., Vollmer et al. 2001), i.e. a 

Virgo cluster has the advantage that it is spiral-rich and j^^g^y eccentric orbit will lead the galaxy at a high veloc- 

dynamically young, i.e. there are strong ongomg environ- ¿^y ¿^^^ ^j^^^^^ ^^^^^ ^^^^^ intracluster médium 

mental mteractions. Chung et al. (2007) found seven spiral ¿^^^^^^ pressure will be very strong. Recently, 

Vollmer (2009) established a model-based ram pressure 
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and the timestep of the best-fit model has to be consistcnt 
with the galaxy's projected position and radial velocity 
within the Virgo cluster. Moreover, Vollmer (2009) con- 
cluded that the linear orbital segments derived from the 
dynamical models together with the intracluster médium 
density distribution derived from X-ray observations are 
consistent with the dynamical simulations of 6 Virgo clus- 
ter spirals. The deduced closest approach of the galaxies 
to the cluster center (M 87) is between 200 and 600 kpc. 

Another consequence of ongoing ram pressure can be 
found in the distribution of polarized radio continuum 
emission. Polarized radio continuum emission is sensitive 
to gas compression and shear motions. Asymmetric ridges 
of polarized radio continuum emission are frequently ob- 
served in Virgo spiral galaxies with disturbed or truncated 
Hi disks (Wezgowiec et al. 2007, Vollmer et al. 2007). 
MHD modelling based on the velocity fields of the dy- 
namical model can corrobórate a ram pressure stripping 
scenario and/or discrimínate between different interac- 
tion scenarios (NGC 4654, Sóida et al. 2006; NGC 4522, 
Vollmer et al. 2006; NGC 4501, Vollmer et al. 2008). 

Despite these clear signs of ram pressure stripping on 
the atomic gas content, only a handful of gas tails asso- 
ciated to cluster spiral galaxies have been detected in X- 
rays: in ACO 2125 (C153 at z = 0.253, Wang et al. 2004), 
in ACO 1367 (UGC 6697 at z = 0.022, Sun & Vikhlinin 
2005), in the Coma cluster (NGC 4848, Finoguenov et al. 
2004), and in ACO 3627 (ESO 137-001, Sun et al. 2006). 
During a ram pressure stripping event it is expected that 
gas is heated to X-ray temperature by bow shocks, heat 
conduction, and the mixing of the stripped ISM into the 
hot intracluster médium (e.g., Stevens et al. 1999, Schulz 
& Struck 2001, Roediger et al. 2006, Roediger & Brüggen 
2008, Tonnesen & Bryan 2009). So it is rather surpris- 
ing and not clear why the X-ray detection rate is so low. 
In particular, there is no detection of the leading bow 
shock of a cluster spiral galaxy in the literature to our 
knowledge (Rasmussen et al. 2006 detected a shock-like 
feature in NGC 2276 which is part of the NGC 2300 group 
of galaxies). Since galaxies move with velocities between 
1000 and 2000 kms~^ and the sound speed in the intra- 
cluster médium of the Virgo cluster is ~ 500-700 kms~^, 
a bow shock with an associated Mach cone is expected. 

In this article we use XMM-Newton data to search for 
X-ray signatures of 3 ram pressure stripped spiral galax- 
ies belonging to the ram pressure stripping sequence from 
Vollmer (2009). We made deep XMM-Newton spectro- 
scopic observations of NGC 4569 and use archival XMM- 
Newton data for NGC 4388, and NGC 4501. They are de- 
scribed in Sect.[2]and the results are presented in Sect.[21 
We find extraplanar diffuse gas in all galaxies. Simple 
Mach cone geometries are fitted to the X-ray images and 
the thermal pressure is compared to ram and cosmic ray 
pressure in Sect. |31 We give our conclusions in Sect. [SI 



2. Observations 

We made deep XMM-Newton observations of NGC 4569 
(proposal ID: 2006501, PI: M. Ehle). In addition, archival 
XMM-Newton observation of NGC 4388 and NGC 4501 
are used. The main astronomical parameters of the stud- 
ied galaxies are presented in Table [1] and the properties 
of their X-ray observations are presented in Table [51 The 
data was processed using the SAS 9.0 package (Gabriel et 
al. I2004P with standard reduction procedures. Following 
the routine of tasks epchain and emchain calibrated event 
lists for the pn camera (Striider et al. , 2001) and for the two 
MOS cameras (Turner et al. I200ip were obtained for each 
galaxy. Next, the event lists were carefully filtered for bad 
CCD pixels and periods of intense radiation of high energy 
background. The filtered event lists were used to produce 
images, background images, exposure maps (without and 
with vignetting correction), masked for an acceptable de- 
tector área using the images scrip10. All images and maps 
were produced in the band of 0.2 - 1 keV in order to be sen- 
sitive for the softest X-ray emission due to hot gas. Next, 
the final images were combined using the data from all 
cameras. The resulting images were adaptively smoothed 
with a Gaussian beam of 10" HPBW. To get a better sig- 
nal to noise ratio and thus better sensitivity for extended 
structures, which allowed to bring up the large-scale struc- 
ture, the images were again smoothed, this time with a 
Gaussian beam of 30" or 1' using the AIPS package. The 
gaps between instrument detectors were not interpolated, 
but filled by combining scaled MOS and pn maps. 

For the spectral analysis, the images were searched for 
point-like sources using the routine task edetectjzhain. 
Next, after excluding circular áreas around the detected 
sources from the relevant región área a spectrum was ac- 
quired. Sizes of circular áreas were determined by visual 
inspection of an image of each point source in the pn event 
list. Blank sky event lists (see Cárter & Read I2007|) were 
used for creating background spectrum corresponding to 
a proper región. For each spectrum response matrices and 
effective área files were produced. For the latter, detector 
maps needed for extended emission analysis were also cre- 
ated. Finally, the spectra were then fitted using XSPEC 11 
(Arnaud 1996). The thermal plasma MEKAL model is 
based on calculations by Mewe and Kaastra (Mewe et 
al íTÜSSl Kaastra [T992| . 

3. Results 

In the following we present the X-ray distribution in the 
band of 0.2 - 1 keV and X-ray spectra of NGC 4569, 
NGC 4388, and NGC 4501 from regions presented in 
Appendix A. The unabsorbed X-ray fiuxes, X-ray lumi- 
nosities, gas densities and temperatures are presented in 
Table [3l The volumes of emitting regions were estimated 
in the following way (see Appendix A) : (i) the galactic disk 
regions are approximated by cylinders with radii set by the 

^ http:/ /xmm.esac.esa.int/external/xmm_science/ 
gallery /utils / images .shtml 
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Table 1. Basic astronomical propcrties of studied galaxies 



Ñame 


Morph. 


Optical position^ 


Incl.^' 


Pos. 


Proj. dist. 






type" 


Q2000 ¿2000 


n 


ang.^n 


to Vir A [°] 


[10^° cm-2] 


NGC 4388 


Sb 


12''25™46':8 +12°39'44" 


82 


91 


1.3 


2.58 


NGC 4501 


Sb 


12''31™59:3 +14°25'14" 


60 


138 


2.0 


2.62 


NGC 4569 


SABa 


12'^36'^5(r.l +13°09'46" 


66 


23 


1.7 


2.82 



taken from HYPERLEDA datábase - |http: / /leda.univ-lyonl.fr| - see Paturel et al. PKÜ]) . 
weighted average valué after LAB Survey of Galactic Hl, see Kalberla et al. ()2005[l . 

Table 2. Parameters of X-ray observations of studied galaxies 



Ñame 


Obs. ID 


Obs. date 


Exp. time'' 


pn filter 


pn obs. 

ñ 


MOS filter 


MOS obs. 

C) 


NGC 4388 


0110930701 


2002-12-12 


11.9 (6.9) 


Thin 


EF 


Médium 


FF 


NGC 4501 


0112550801 


2001-12-04 


14 (2.9) 


Médium 


EF 


Thin 


FF 


NGC 4569 


0200650101 


2004-12-13/14 


66 (49) 


Thin 


EF 


Médium 


FF 



^ Total time in ksec with clean time for pn camera in brackets. 
Observing mode: FF - FuU Erame, EF - Extended FuU Frame. 



Table 4. Estimated volumes and gas masses. 

Galaxy Región volume (mI^^V'^T (M^lrf^ 

[pe"] [Me] [Me] 

NGC disk TTTeTTÍ e.Gttle+Oñ 7.8+4>-|-06 

4569 halo 6.14e-H2 9.8Í|?e-K07 1.0Íg;^e+08 

NGC nucleus LlOe+To a9TJfe+06 4.7+2.86+06 

4388 disk 4.52e-fll 4.6Í2:7e+07 3.7t2:}e+07 

tail 1.38e-H3 2.1í};^e-K08 2.8t|ie-|-08 

NGC disk 7T0e+TÍ - 5.4+i^;^e-H07 

4501 tail 1.91e-H3 - 2.5tile+08 

rj is the volume fiUing factor. 



extent of X-ray emission and a height of 1 kpc, (ii) the halo 
emission of NGC 4569 is assumed to be spherical with a 
radius of 11.4 kpc, (iii) the tail región of NGC 4388 is ap- 
proximated by a cylinder with a length of 50.8 kpc and a 
radius of 9.3 kpc, and (iv) the tail región of NGC 4501 is 
approximated by a cylinder with a length of 30 kpc and 
a radius of 14 kpc. The model types and reduced ^-re 
shown in Table [5] The volumes are uncertain by roughly 
a factor of 2. The typical uncertainty of the density is 
30 — 50%, that of the masses is ~ 0.2 dex. The volumes 
and gas masses of the different X-ray emitting components 
are presented in Table Hl 

3.1. NGC 4569 

The difFuse X-ray emission associated to NGC 4569 is ex- 
tended over an área of approximately 20' x 20' (100 x 
100 kpc; Fig. [T]). The overall shape of the X-ray emission 
is that of a triangle with the tip pointing to the northeast. 



The galactic disk has a 5 times lower X-ray ñux than the 
giant difFuse halo which encompasses also a dwarf satellite 
galaxy IC 3583. The large radio lobes observed at 20 cm 
(Chyzy et al. 2006) are located entirely within the dif- 
fuse X-ray halo. There is enhanced X-ray emission asso- 
ciated to the Southern half of the western radio lobe (at 
RA=12h36m30s DEC==13°09'). 

X-ray spectra (Fig. [2]) have been taken from the disk 
and the inner halo (2.5' x 2.5' or 12 x 12 kpc) excluding 
emission from the galactic disk. They are fitted by a two 
temperature thermal plasma MEKAL model with a hot 
and a cold component. An additional powerlaw component 
is added to account for undetected point sources. Both, 
disk and halo show a cold (~ 0.15 keV) and a hot com- 
ponent (~ 0.5 keV) and densities oí n ^ 10"'^ 77"°'^ cm""^ 
(Table [3]), where rj is the volume filling factor. 

3.2. NGC 4388 

The galactic disk of NGC 4388 is detected in X-rays. It is 
more extended to the north than to the south. In addition, 
an extended tail is detected towards the northeast of the 
galactic disk. This X-ray tail is associated with the Ha 
plume reported by Yoshida et al. (|2002p and the beginning 
of the huge Hi tail discovered by Oosterloo & van Gorkom 
(2005) (Fig. [3]). At distances > 5' from the galactic disk 
the X-ray emission starts to be dominated by the bright 
halo of M 86. 

We performed a spectral analysis of the nucleus, the 
disk and nucleus, and the part of the tail cióse to the 
galactic disk. The spectra (Fig. S]) are fitted with two sin- 
gle temperature thermal plasma MEKAL models with an 
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Table 3. X-ray fluxes, derived densities, and temperatures from the X-ray spectra. 



Galaxy 


Región 


hot gas flux^ 


Lx 


(ncoid v"-'r 


^^cold 










[10"^'' ergcm-^s-i] 


[lO^^'ergs-i] 


[10"^cm~^] 


[keV] 


[10"^cm-^] 


[keV] 


NGC 
4569 


disk 
halo 


2.26ÍÍ:- 


0.780 
3.936 


Q 7+0.3 


14+"1" 
17+002 

"•^'-0.03 


7+°-2 


0.46±0.08 
0.49±0.05 


NGC 


nucleus 


16.84t-.- 


5.819 


r)rT Q + Í3.Y 

^'•«-10. 6 




18.81^;^ 


O.66I0.04 


4388 


disk 


10 44+'*-'^'^ 


3.608 




0.21±0.03 


o 7+1.0 
•J- ' -1.0 


0.79±0.11 




tail 


21.51_g g2 


7.433 




0.87+°;}^ 


9+^ '^ 

"•^-0.9 


2.26í^i« 


NGC 
4501 


disk 
tail 


16.44t^:S¿ 


5.681 
4.945 






n . + 1.U 
■^•^-1.0 

fi+°-2 

U.O_Q 2 


o.38Í^:- 



Total unabsorbed flux from the thermal emission between 0.3 and 12 keV. 



77 is the volume fillirig factor. 



Table 5. Model type and reduced Xri 



Galaxy 


model type 


Xred 


NGC 4569 disk 


wabs(mekal+mekal+powerlaw) 


0.96 


NGC 4569 halo 


wabs(mekal+mekal+powerlaw) 


1.23 


NGC 4388 nucleus 


wabs(mekal+mekal+powerlaw+wabs(powerlaw+gauss)) 


1.06 


NGC 4388 disk+nucleus 


wabs(mekal+mekal+powerlaw+wabs(powerlaw+gauss)) 


1.07 


NGC 4388 tail 


wabs(mekal+mekal+powerlaw) 


1.11 


NGC 4501 disk 


wabs(mekal+powerlaw) 


1.39 


NGC 4501 tail 


wabs(mekal+powerlaw) 


0.46 



additional powerlaw component to account for undetected 
point sources. In the fit of the disk and the nucleus, a con- 
tribution from the AGN was included by adding to the 
model an absorbed powerlaw component and the 6.4 keV 
iron line (see Table [5]) . The temperatures and densities 
are given in Table |31 The nucleus has a gas density of 
~ 30?7~°'^ cm~'^. The cold and hot components of the 
disk have temperatures of ~ 0.2 keV and ~ 0.8 keV, re- 
spectively. The exceptionally high temperature of the hot 
component in the disk is probably due to the nuclear out- 
flow (Veilleux et al. 1999, Yoshida et al. 2004). In the tail 
the temperature of the cold component is comparable to 
the hot component of the disk. The hot component cor- 
responds to the intracluster médium. About 60 % of the 
emission comes from the hotter spectral component. 



3.3. NGC 4501 

NGC 4501 shows a moderately extended hot gas halo 
(Fig. [5]). The northeastern side is more extended than the 
southwestern side. X-ray emission is found up to ^ 3' 
15 kpc) from the galactic disk north from the galaxy 
center. This X-ray plume is cióse to the Ha-bright spiral 
arm at the outer northwestern disk which has no counter- 
part in the southeast. 

Due to the limited number of counts, the X-ray spectra 
of the disk and the tail (Fig. |6]) has been fitted by a single 
hot component of ^ 0.4 keV and ~ 0.7 keV, respectively 
(Table H. 



4. Discussion 

To put the observed X-ray emission and spectra of the 3 
galaxies into context, we first give a brief summary on the 
evolutionary stage of each galaxy. 

4.1. NGC 4569 

The spiral galaxy NGC 4569 is a very peculiar member of 
the Virgo cluster. Its very large size {D25 = 9.5', or 47 kpc) 
compared to other cluster members, and negative radial 
velocity (ui- — —235 kms~^ with respect to the cluster 
mean velocity of ~ 1100 kms~^) have lead to some doubt 
about its cluster membership, despite its cióse projected 
distance to the cluster center {d — 1.7° = 505 kpc). Van 
den Bergh (1976) classified this very red galaxy as anemic 
due to its low arm-interarm contrast on optical images. 
Hi line synthesis observations (Warmels 1988; Cayatte et 
al. 1990) showed that it has lost more than 90% of its 
Hi gas and that its Hi distribution is heavily truncated. 
NGC 4569 is thus an exceptional galaxy, it is one of the 
largest and most Hi deficient galaxies in the Virgo cluster 
(Solanes et al. 2001). 

More recent Hi observations showed a low surface den- 
sity arm in the west of the galaxy, whose velocity field is 
distinct from that of the overall disk rotation. The ob- 
served gas distribution, velocity field, and velocity disper- 
sión are consistent with a major stripping event that took 
place about 300 Myr ago. This post-stripping scenario can 
reproduce the main observed characteristics of NGC 4569 
(Vollmer et al. 2004). The analysis of the multiwavelength 
photometry of the galactic disk including GALEX UV ob- 
servations (Boselli et al. 2006) confirmed this scenario. 
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NGC 4Seg XMM Newton M0S1+M032+pn 0.2 1 keV on DSS blue 



NGC456Í) disk 




13 37 30 



00 3645 30 1S 

RIGHT ASCENSION (jaOOO) 

NGC 4563 VLA 1 .4 GHz TP on XMM-Hewton MOSI 4*<OSa4^pn OS - t keV 



LU 
O 




00 36 45 30 15 

RIGHT ASCENSION (J2a00) 

Fig. 1. X-ray emission (0.2-1 keV) of NGC 4569 together 
with an optical B band (upper panel) and the radio total 
power map at 1.4 GHz froni Chyzy et al. (¡20061 lower 
panel). The resolution of both maps is 1'. The beam size 
is shown in the bottom left cerner of the figure. The first 
greyscale in the lower panel corresponds to 3 times the 
rms noise level. The radio contour levéis are 3, 5, 8, 16, 
25, 40, 80, 120x0.24 mJy. 



Tschóke et al. (2001) found a diffuse soft X-ray component 
extended to the west coinciding with a giant Ha structure. 
They concluded that the X-ray gas traces a large scale out- 
ñow from accumulating supernova explosiona and stellar 
winds in the galactic center. 

Surprisingly, large radio continuum lobes were found in 
NGC 4569 by Chyzy et al. (2006). These lobes extend up 
to 24 kpc from the galactic disk and are the first observed 
in a cluster spiral galaxy. Since the radio lobes are rather 
symmetric, they appear to resist ram pressure. The most 




chamicl energy (keV ) 



NGC4569 halo 





t 








^ 



channel energy (keV) 



Fig. 2. X-ray spectrum of NGC 4569 's disk (upper panel) 
and inner halo excluding emission from the galactic disk 
(lower panel) . 



plausible explanation for the origin of the radio lobes is 
a galactic superwind. To resist ram pressure, the pressure 
within the región of the radio lobes has to be higher than 
ram pressure. Assuming the mínimum energy condition, 
the total pressure of cosmic rays and magnetic fields in 
the lobes is p ^ 10~^^ dyncm"^. With a mean outflow 
velocity of Wwind — 700 kms~^ the mean density is n ^ 
3 X 10~* cm~^ (Chyzy et al. 2006). A sufficiently strong 
superwind-type flow can easily overeóme the gravitational 
potential of the galaxy and resist ram pressure. The fact 
that the radio lobes are entirely engulfed by the diffuse 
X-ray halo of NGC 4569 lends support to this scenario. 

The thermal pressure from the halo X-ray spectrum 
(Fig. [2]) can be compared to the pressure of the cosmic 
rays and magnetic fields from Chyzy et al. (2006). We 
find a pressure of p^^''^ ^ 2 x lO^^'^yy^"'^ dyncm^^ for 
the cold component (T = 2 x 10^ K) and pj^'^* - 5.5 x 
10~^'^77~'^'^ dyncm^^ for the hot component (T — 5.7 x 
10® K), where rj is the volume fiUing factor. The thermal 
pressure of the hot phasc is thus comparable, but a factor 
of 2 lower than that of the cosmic rays and magnetic fields. 
Thus the superwind might be driven by cosmic rays. 
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HOC 4388 XMM-Newlon MOSI+MOSa+pn O.g- 1 kBV on DSS Mué 

laso 



NGC 4388 nuclcus 



g 
1- 




mam ascenision <J2üop} 

JWaCáaaSH IbH XMIa^Navuten MC>S1»MOSai-pn 0.g - 1 hétf 



36 



1300 




30 15 00 

nGHV ASCENSION (JEOCO) 



Fig. 3. X-ray emission (0.2-1 keV) of NGC 4388 together 
with an optical B band image (upper panel) and the Hi 
map from Oosterloo & van Gorkom (2005; lower panel). 
The X-ray map was cut at the position, where contribution 
from the M 86 halo starts to domínate. The first greyscale 
corresponda to 7 times the rms noise level. 



The triangular shape of NGC 4569's diffuse X-ray 

halo reminds that of a bow shock. From the dynami- 
cal model of Vollmer et al. (2005) and Vollmer (2009) 
we know the direction of the galaxies motion within the 



— -3 p 







1 J 


- 








j 1 1 1 : 







1 2 

channel energy (keV) 



I o 




3 o 

a 



tfti 



0.5 

NGC4388 taU 



1 2 5 

channel energy (keV) 




0.5 1 

channel energy (keV) 



Fig. 4. The X-ray spectra of NGC 4388's nucleus (upper 
panel), disk and nucleus (middle panel), and tail (lower 
panel). 



intracluster médium. With an estimated galaxy velocity 

of Ugai = 1500 kms"-'^ and a sound speed of c = 550- 
700 kms^^ (depending on the adiabatic index) in the in- 
tracluster médium the Mach number is M 2.1-2.7 and 
the angle of the Mach cone (half of the opcning angle) is 
a = sin~^(-¡g) = 22-28°. On the other hand, since the very 
tenuous (n 10~^ cm~^) intracluster médium is fully ion- 
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NGC 4501 XMM-Newton M0SUM0S2+pn 0.2 ■ 1 keV on DSS blue 



NGC4501 disk 



O 
< 



14 32 



30 




12 32 20 15 
NGC 4501 XMM-Newton M0S1+M0S2+pn 0.2 - 1 keV on H alpha 



O 
I- 
< 



O 
UJ 
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26 



channel energy (keV) 



12 32 25 



15 



10 05 00 31 55 50 
RIGHT ASCENSION (J2000) 



45 



40 



Fig. 5. X-ray emission (0.2-1 keV) of NGC 4501 overlaid 
onto an optical B band image (upper panel) and the Ha 
image from Goldmine (lower panel). The map resolution 
is 30". The beam size is shown in the bottom right córner 
of the figure. The X-ray contours are 3, 5, 8, 16, 25 times 
the rms noise level. 



ized and magnetized, the sound speed has to be replaced 
by the magnetosonic velocity v¡ns = \/ (c^ -I- v'jj, where 
the Alfvén speed is v¡¡. = , ^ Typical magnetic field 

w47rnmp 

strengths in galaxy clusters are of the order of a few /íG 
(Ferrari et al. 2008). 

We overlaid a Mach cone on the X-ray emission dis- 
tribution of NGC 4569 in Fig. [71 The direction is based 



Fig. 6. The X-ray spectra of NGC 450rs disk (upper 
panel) and tail (lower panel). 



on the dynamical model whereas the opening angle is fit- 
ted by eye to the triangular form of the X-ray distribu- 
tion. In this way we find an apparent opening angle of 
75° or an apparent Mach cone angle of a = 37.5°. The 
true Mach cone angle depends on the projection under 
which the Mach cone is observed, i.e. the angle /3 between 
the galaxy's normalized 3D velocity vector and the plañe 
of the sky. Given the normalized 3D model velocity vec- 
tor (-0.64,0.53,0.56) (Vollmer 2009), this angle is /? = 34° 
and the true Mach cone angle is 30°. For (3 = 20°/40° 
the true Mach cone angles are 36° and 28°, respectively. 
NGC 4569's high radial velocity with respect to the Virgo 
cluster mean favors higher valúes of /3. The true Mach cone 
angle is thus consistent, but might be somewhat higher 
than our estimates based on the sound speed of the intr- 
acluster médium. We thus suggest that the magnetic field 
plays a role and that the magnetosonic velocity is relevant 
for the determination of the Mach number. Based on true 
Mach cone angles between 30° and 32°, a sound speed 
of 550-700 kms"^, a galaxy velocity of 1500 kms~^, and 
assuming an intracluster density of n = 10~^ cm~^, we 
derive Alfvén speeds between 300 kms^^ and 600 kms— 1 
and a magnetic field strength of ^ 1-3 fiG at the loca- 
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tion of NGC 4569. The derived range of Alfvén speeds is 
thus comparable to that of the sound speed. The derived 
magnetic field strength is consistent with typical magnetic 
field strengths in galaxy clusters (Ferrari et al. 2008). 




12^ 38^ 00^37"^ 30°' 00^ 36^ 30^ 00^ 35^ 30^ 
Right Ascensión 
Center: R.A. 12 36 49.01 Dec +13 11 33.1 

Fig. 7. X-ray emission (0.2-1 keV) of NGC 4569 together 
with the Mach cone whose direction is based on the dy- 
namical model and whose opening angle is adjusted to the 
X-ray distribution. 



4.2. NGC 4388 

NGC 4388 's projected location within the Virgo cluster is 
cióse to M 87 (1.3° = 385 kpc). Yoshida et al. (|2DD4| 
discovered an Ha plume which extends up to 35 kpc 
northeast from the galaxy center. Based on a Effelsberg 
lOOm telescope Hi detection of the Ha plume, VoUmer & 
Huchtmeier (2003) made dynamical models and concluded 
that ram pressure stripping can account for the observed 
extraplanar Ha and Hi emission. Oosterloo & van Gorkom 
(2005) observed NGC 4388 with the WSRT and discov- 
ered an Hi tail which extends more than 100 kpc to the 
north. 

We clearly detect X-ray emission associated to the Hi 
and Ha tail of NGC 4388. We can only attribute the ob- 
served X-ray emission to the first ^ 5' (^ 25 kpc) of the 
Hi tail which is about 5 times longer. This is because the 
projected Hl distribution approaches M 86 and the halo 
of M 86 dominates the X-ray emission in this región. 

As for NGC 4569 we adjusted a Mach cone to the X- 
ray and Hl emission distribution with the direction being 
based on the dynamical model (Vollmer & Huchtmeier 
2003, Vollmer 2009). The apparent Mach cone angle is 
a ^ 26° (Fig. [5|). With a normalized 3D model velocity 



vector (0.09,-0.75,0.66) (Vollmer 2009) the true Mach cone 
angle is 20°. This corresponds to a Mach number of 2.9 
and a galaxy velocity of ^ 1600-2000 km s~^ (with a sound 
speed of 550-700 kms^^), which is consistent with the 
estimated galaxy velocity based on the dynamical model 
(1700 kms"\ Vollmer 2009). To obtain an Alfvénic Mach 
number of 2.9 with Vgai = 1700 kms~^, a magnetic field 
strength of ~ 3 ¡iG is needed with an intracluster density 
of n = 10~^ cm~^. 




12^ 26^ 20^^ 10^ 00^ 25^ 50^ 40^^ 

Right Ascensión 
Center; R.A. 12 25 57.33 Dec +12 42 33.9 



Fig. 8. X-ray emission (0.2-1 keV) of NGC 4388 together 
with the Mach cone whose direction is based on the dy- 
namical model and whose opening angle is adjusted to the 
X-ray distribution. 

Within the X-ray tail we measure a temperature of 
0.87 keV (10"^ K) and a density of n - r]-°-^ 10'^ cir-^. 
The thermal pressure is thus pth = 1-4 x 10~^^ dyncm"^, 
whereas ram pressure is about 4-8 X 10^^^ dyn cm ^. As 
in the halo of NGC 4569, the thermal pressure based on 
the X-ray observations is significantly lower than the esti- 
mated ram pressure acting on NGC 4388. 

We made an attempt to calcúlate the total mass of 
the hot X-ray gas within the whole Hi tail. Unfortunately, 
this is not directly possible due to the dominating X-ray 
halo of M 86. Several assumptions have to be made: first, 
a constant surface-brightness along the whole tail is as- 
sumed. Next, the flux from the área of detectable X-ray 
tail is rescaled to the total volume of the Hl tail. Since the 
Hi tail is ~ 3 times longer than the X-ray tail, this leads 
to a 3 times larger volume. 

Givcn the temperature from the spectral fit, the den- 
sity of the hot gas is n ^ 77"°-^ 10"'^ cm~'^. A total gas 
mass of ~ 2 X 10^ Mq is found in the part of the X- 
ray tail which is not dominated by M 86 's halo (south 
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to the declination 12 48 00). Since our volume estímate 
is on the high side and since we expect a volume filling 
factor smaller than 1, a more realistic mass estímate of 
the X-ray tail is ^ 10* Mq. Within a 3 times longer tail, 
there might be a few 10® Mq of hot gas associated with 
the entire Hi tail of NGC 4388. This total mass of hot 
gas is comparable to the total mass of atomic hydrogen in 
the tail (3.4 x 10*^ M©; Oosterloo & van Gorkom 2005). 
The total gas mass of the tail (~ 6 x 10® Mq) is much 
lower than the expected stripped gas mass based on the 
Hl deficiency (2x10*^1^0; VoUmer & Huchtmeier 2007). 

Possible reasons for this discrepancy are strong ISM- 
ICM mixing with a strong decrease of the gas density or 
star formation within the gas tail. Based on the recon- 
structed star formation history from deep optical spec- 
tra Pappalardo et al. (2010) found that the stripped gas 
left the galactic disk ^ 200 Myr ago. The gas which has 
been stripped from the location where the optical spectra 
were taken is presumably now in middle of the Hi tail. 
The gas within the tail thus traveled with a velocity in 
the sky plañe of ^ 300 kms^^. With a radial velocity of 
^ 300 kms^^ the total outflow velocity of the stripped gas 
is ~ 400-450 kms^^. This 3D configuration is consistent 
with the VoUmer & Huchtmeier (2003) model. The strip- 
ping timescale is reasonable in a scenario where NGC 4388 
has been stripped by a recent passage through the Virgo 
cluster core, cióse to M 87 (VoUmer 2009). Therefore, 
we think that the tail traces the entire stripping process. 
Given the sniall number of Hii regions within the Hi tail 
(Yoshida et al. 2004), gas loss via star formation seems 
not to be significant. 

VoUmer & Huchtmeier (2007) suggested an alternative 
explanation for the missing stripped gas: spiral galaxies 
entering a cluster may rapidly lose the atomic gas which 
is located beyond the optical radius. This gas is (i) less 
bound than the gas in the inner disk and (ii) no longer 
stirred by strong interstellar turbulence and can easily be 
heated and evaporated. The stripping of the outer gas 
disk might happen before the galaxy enters the cluster 
core. In this case the galaxy already shows an initial Hi 
deficiency of 0.4 before the stripping event. If such an ini- 
tial Hi deficiency is assumed, the sum of the observed Hi 
and extrapolated X-ray gas mass might account for the 
whole stripped gas. Alternatively, a fraction of the gas 
mass might be in the form of tenuous ionized gas with 
temperatures below 10^ K. 

4.3. NGC 4501 



that ram pressure stripping can account for the main ob- 
served characteristics of NGC 4501 (VoUmer et al. 2008). 
The galaxy is stripped nearly edge-on, is heading South- 
west, and is 200-300 Myr before peak ram pressure, i.e. its 
closest approach to M 87. 

A possible Mach cone is adjusted to the X-ray emission 
distribution in Fig. [9] The apparent Mach cone angle is 
^ 30°. With a normalized 3D model velocity vector (0.48,- 
0.44,0.76) (VoUmer 2009) the true Mach cone angle is 20° 
implying a Mach number of M ~ 2.9. Assuming a lower 
angle between the galaxy's normalized 3D velocity vector 
and the plañe of the sky (/? = 37° instead of 49° yields a 
true Mach cone angle of 25° and a Mach number of M ~ 
2.4. If the sound speed is relevant for the Mach number 
the galaxy velocity is Wgai ^ 1300-2000 kms^^. This is 
consistent with the galaxy velocity based on the dynamical 
model (tigai = 1500 kms~^). If the Alfvén velocity exceeds 
the sound speed, a magnetic field strength of ^ 2-3 /íG is 
needed assuming an intracluster médium density of n = 
10--* cm-3. 




12" 32^ 30^ 20^ 

Center: R.A. 12 31 57,78 



O 00 31 50' 

Right Ascensión 



40 30 
Dec +14 26 50.4 



Fig. 9. X-ray emission (0.2-1 keV) of NGC 4501 together 
with the Mach cone whose direction is based on the dy- 
namical model and whose opening angle is adjusted to the 
X-ray distribution. 



Polarized radio emission studies showed a strong compres- 
sion of the magnetic field at the southwestern edge of the 
galactic disk (Wezgowiec et al. 120071 VoUmer et al. l2007p . 
NGC 4501 's Hi disk is sharply truncated to the South- 
west, well within the stellar disk. On the opposite side, 
there is low surface-density gas, which is more extended 
than the main Hi disk. A detailed comparison of obser- 
vations with a saniple of dynamical simulations, showed 



Our X-ray data do not provide enough counts to ob- 
tain a high signal-to-noise spectrum. Therefore, we fitted 
a single temperature model to the X-ray spectrum of the 
tail and found the temperature to be ^ 0.7 keV. This 
temperature is comparable to the temperatures of the hot 
extraplanar gas in NGC 4569 (0.49 keV) and NGC 4388 
(0.87 keV). The X-ray emission beyond the Hi disk might 
be due to stripped atomic gas mixed with the intracluster 
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médium or a galactic wind due to star formation in the 
prominent northwcstern spiral arm (Fig. [5]). 

4.4. Hot gas within the Mach cones 

For NGC 4569 we have the lucky situation that the Mach 
cone is entirely fiUed by a galactic outñow. A galactic 
outñow encomitering the ICM ñow within the Mach cone 
should produce an observable contact discontinuity (see, 
e.g., Heinz et al. 2003, Kraft et al. 2011) in addition 
to the bow shock. According to Kraft et al. (2011) and 
Vikhlinin et al. (2001) the separation between the two 
features is ^ 0.15 times the distance between the galaxy 
center and the contact discontinuity for large Mach num- 
bers. Assuming this distance to be 25 kpc, the projected 
distance between the bow shock and a contact discontinu- 
ity is about 45". This is cióse to the image resolution of 
1'. The fact that we do not observe two distinct features 
might be caused by the absence of a contact discontinuity, 
because the ICM ñow within the Mach cone is turbulent, 
or might be due to resolution/sensitivity and/or projec- 
tion effects. 

On the other hand, the Mach cone geometries of 
NGC 4388 and NGC 4501 are tentative. Ah extraplanar 
X-ray emitting hot gas is found within the Mach cones of 
the three galaxies. 

Whereas it is most likely that the diffuse X-ray halo of 
NGC 4569 (Fig.[T|) is due to a galactic superwind originat- 
ing in a nuclear starburst, the extraplanar X-ray emission 
of NGC 4388 is most likely due to the mixing of stripped 
ISM with the hot intracluster médium. Typical temper- 
atures of the cold and hot phase in the star formation 
driven outñow of NGC 253 are ~ 0.1 keV and ~ 0.3 keV, 
respectively (e.g. Bauer et al. 2008). The appearance of 
a "multi-temperature" halo is a natural consequence of 
the plasma driven out of collisional ionization equilibrium 
(CIE) towards a non-equilibrium ionisation (NEI) state in 
a galactic outñow. The temperature of the hottest com- 
ponent in galactic outñows of massive spiral galaxies is 
about 0.4 keV (Martin 1999). 

The gas in NGC 4388's and NGC 4501 's tails is sig- 
nificantly hotter than that of typical galactic outñows of 
massive spiral galaxies. On the other hand, the X-ray tail 
of ESO 137-001 in ACO 3627 has a temperature of 0.7 keV 
(Sun et al. 2006). The temperature of the X-ray tail of 
ESO 137-002 might be as high as 2 keV (Sun et al. 2010). 
Both galaxies evolve in a galaxy cluster whose intracluster 
médium is 3 times hotter than that of the Virgo cluster 
(2.4 keV, Boringer et al. 1994). We suggest that the higher 
gas temperatures observed in the X-ray tails of NGC 4388 
(0.9 keV) and NGC 4501 (0.7 keV) are due to the mixing 
of the stripped ISM into the hot intracluster médium of 
the Virgo cluster. 

5. Conclusions 

XMM-Newton observations of three Virgo cluster spirals 
are presented. These galaxies are part of the ram pres- 



sure stripping time sequence of Vollmer (2009): NGC 4501 
is still approaching the cluster center, NGC 4388 and 
NGC 4569 have experienced peak ram pressure 150 Myr 
ago and ~ 300 Myr ago. We ñnd extraplanar diffuse X-ray 
emission in all galaxies: 

— In NGC 4569 it is associated with the galactic out- 
ñow observed in Ha and radio continuum (Chyzy et 
al. 2006). The diffuse X-ray emission has a triangular 
shape. 

— In NGC 4388 the X-ray emission follows the beginning 
of the Hi tail discovered by Oosterloo & van Gorkom 
(2003). The halo of M 86 dominates the X-ray emission 
at larger distances from the galaxy. 

— In NGC 4501 it is located opposite to the side of the 
galactic disk where the gas is compressed by ram pres- 
sure (Vollmer et al. 2008). 

Based on the 3D velocity vectors from dynamical mod- 
elling (Vollmer 2009) a simple Mach cone can be ñtted to 
the triangular shape of NGC 4569 's diffuse X-ray emis- 
sion. To our knowledge this is the ñrst direct evidence for 
a Mach cone in a cluster spiral galaxy. We can only ob- 
serve it, because it is ñlled with hot gas from a galactic 
superwind. Assuming that all extraplanar diffuse X-ray 
emission has to be located inside the Mach cone, we also 
fit Mach cones to NGC 4388 and NGC 4501. 

For NGC 4569 it is hard to reconcile the derived Mach 
cone opening angle with a Mach number based on the 
sound speed alone. Instead, a Mach number involving the 
Alfvénic speed seems to be more appropriate, yielding a 
magnetic ñeld strength of ~ 3-6 /íG for a intracluster 
médium density of n ^ 10^^ cm~^. The shock is thus 
caused by fast magnetosonic waves. 

Gas densities and temperatures are derived from the 
X-ray spectra. In the outñow of NGC 4569 the thermal 
pressure is a factor of 2 smaller than the cosmic ray gas 
pressure. The outñow might thus be driven by cosmic rays. 
The thermal pressure in NGC 4388's X-ray tail is signiñ- 
cantly smaller than the ram pressure which NGC 4388 is 
experiencing. Whereas the temperature of the hot compo- 
nent of NGC 4569's X-ray halo (0.5 keV) is at the high 
end but typical for a galactic outñow, the temperature of 
the hot gas tails of NGC 4388 and NGC 4501 are signifi- 
cantly hotter (0.7-0.9 keV). We suggest that these higher 
gas temperatures are due to the mixing of the stripped 
ISM into the hot intracluster médium of the Virgo clus- 
ter. 
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Appendix A: X-ray emitting volumes 
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Fig. A.l. Regions for X-ray spectrum extraction and X- 
ray emitting volume calculations. 



